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ABSTRACT 

, We present a theoretical model for the alternating phase lags associated with 

QPO fundamental and harmonic frequencies observed in some Galactic black-hole 
candidates. We assume that the accretion flow exhibits a transition from an outer 
cool, optically thick accretion disk to a hot, inner advection-dominated accretion flow 
(ADAF), and that the QPOs are related to small-scale oscillations of the accretion 
rate and the transition radius. We present an analytical estimate of the expected 
phase lags at the fundamental and first harmonic frequencies of the QPOs and perform 
detailed time-dependent Monte-Carlo simualtions of the radiation transport in the 
oscillating ADAF / cool disk system. We find that this model is well suited to 
\ reproduce alternating phase lags between the fundamental and the first harmonic. It 

also naturally explains the trend observed in GRS 1915+105 that, as the soft X-ray 
' luminosity increases, the QPO frequency increases and the phase lag associated with 

the QPO fundamental frequency changes sign from positive to negative. The relation 
between the disk temperature and the QPO frequency observed in GRS 1915+105 is 
consistent with a secular instability modulating the disk evaporation at the transition 
q ! radius. 

CO 

a 
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radiative transfer — radiation mechanisms: thermal 

1. Introduction 

The steadily increasing amount of high-quality, high time-resolution X-ray data from Galactic 
X-ray binaries, has stimulated vital interest in the characteristics of the rapid variability of these 
objects. The X-ray emission from both Galactic black-hole candidates (GBHCs) and accreting 
neutron stars is known to vary on a wide range of time scales, sometimes showing quasi-periodic 



oscillations (QPOs) (for reviews see, e.g., |van der Klis 1995] , |Cui 1999a|) . In GBHCs, such QPOs 



have been detected at frequencies ranging from several mHz to <^ 100 Hz. These QPOs are most 
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notable when the sources are in the low-hard or in the rare very high state. Recently, correlations 
between the amplitudes and centroid frequencies of several types of QPOs with the spectral 
characteristics of the source have been found in some objects ( Rutledge et al, 1999| , Markwardt, 



Swank, & Taam 1999, Sobczak et al. 1999, Muno, Morgan & Remillard 1999). A general pattern 
emerging from these analyses is that the frequency of those types of QPOs which do show a 
correlation with spectral properties, seems to increase with both the power-law and the disk 
black-body flux from the source. 

A surprising property of some types of QPOs has recently been found in the 67 mHz QPO of 
GRS 1915+105 (|Cui 1999b!) , the 0.5 - 10 Hz QPO of the same source (|Reig et al. 2000| , |Lin et al. 



2000[) , and in the 0.3 - 3 Hz QPO in XTE J1550-564 (Cui, Zhang, & Chen 2000): While the phase 
lag at the QPO fundamental frequency is negative, the phase lag associated with the first harmonic 
was positive. In the case of the low-frequency QPO of GRS 1915+105, even three harmonics were 
detected, and the phase lags were found to alternate between subsequent harmonics ( |Cui 1999b| ). 
The phase lag associated with the 0.5 - 10 Hz QPO of GRS 1915+105 was found to change sign 
from positive to negative as the QPO frequency increases above ~ 2.5 Hz ( [Reig et al. ^000| , |Lin et| 
20001 ). 



al 



These peculiar patterns are apparently completely counter-intuitive in the light of currently 
discussed models for the hard phase lags in X-ray binaries. Models proposed to explain the hard 
phase lags are either based to the energy-dependent photon escape time in Compton upscattering 
scenarios for the production of hard X-rays (Kazanas, Hua fc Titarchuk 1997, Bottcher & Liang 
1998), or due to intrinsic spectral hardening during X-ray flares, e.g. due to decreasing Compton 
cooling in active regions pushed away from an underlying accretion disk in a patchy-corona model 



(Poutanen & Fabian 1999) or due to density perturbations drifting inward through an ADAF 
toward the event horizon (Bottcher & Liang 1999). These models dealt only with the continuum 
variability and did not consider the effects of QPOs. 

In this paper we will explore the response of a two-phase accretion flow, consisting of an outer, 
cool, optically thick accretion disk and an inner, hot ADAF (Narayan & Yi 1994, Abramowicz et 
al. 1995| , |Chen et al. 1995 ), to a periodically varying soft photon input from the cool disk ( [Liang 
fc Bottcher 200C| ). A two-phase accretion flow with an inner ADAF has been found to produce 
good fits to the photon spectra of, e.g., several Galactic X-ray binaries (e.g., Narayan, McClintock 
& Yi 1996, |Hameury et al. 19971 , |Esin et al. 1998[) , low-luminosity AGN (|Quataert et al. 1999| ) 
and giant elliptical galaxies ( [Fabian fc Rees 1995 , di Matteo fc Fabian 1997 , di Matteo et al.| 
2000 ). In §||, we describe the basic model setup according to this two-phase accretion flow, and 
derive an analytical estimate for the expected phase lags associated with the QPO and the first 
harmonic applicable in some simplified cases. A short description of the Monte-Carlo simulations 
used to solve the time-dependent radiation transport problem follows in §||. In §|| we describe a 
series of simulations designed specifically to explain the peculiar phase lag behavior associated 
with the 0.5 - 10 Hz QPO in GRS 1915+105. We summarize in §§. 
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2. Model setup and analytical estimates 

The basic model setup, as motivated in the introduction, assumes that the inner portion of 
the accretion flow is described by an ADAF, where the density profile is approximately given by 
a free-fall profile, n e (r) oc r -3 / 2 , and the electron temperature is close to the virial temperature 
and thus scales as T e (r) oc r _1 . We emphasize, however, that other hot two-temperature 
inner disk models would give similar results. The ADAF exists out to a transition radius 
-Rtn beyond which the flow is organized in a standard optically thick, geometrically thin 



Shakura-Sunyaev disk (Shakura fc Sunyaev 1972). The transition radius is typically expected to 



be several 100 R s £ R tT £ 10 4 R s ( ponma 199(j |Manmoto et al. 200q , |Meyer et al. 2000|) , where 



R s is the Schwarzschild radius. Outside this transition radius, the disk temperature scales as 
To(r) oc r~ 3 / 4 . With the free-fall density profile given above, the radial Thomson depth of the 
corona may be written as t adaf = 0.71 Mn / {mr\^ 2 ), where Myj is the accretion rate through 
the ADAF in units of 10 17 g s , m is the mass of the central compact object in units of solar 
masses, and rj is the inner edge of the ADAF in units of Schwarzschild radii. 

The basic assumption of our baseline model is that the observed QPOs in the X-ray variability 
of Galactic black-hole candidates are related to small-scale oscillations of the transition radius. 
Such oscillations in r tr may be caused by variations of the accretion rate, but we defer a detailed 
study of possible mechanisms driving the oscillations to a later paper. We thus assume that, as a 
function of time t, the transition radius oscillates as Rtr(t) = r tr (1 + o, T sin[cji]), where a r <C 1 is 
the amplitude of the oscillation and oj = 27t/qpo the QPO frequency. Denoting £ = 1 + a r sin(cji), 
the time-dependent disk flux may then be estimated as Fj)(t) oc [Rt r (t)] 2 [Trf(t)] 4 oc (r tr £) _1 . 
In the following, we are focusing on an analytical description of the X-ray signals at the QPO 
fundamental and first harmonic frequencies. Thus, we will expand all expressions up to the 2nd 
order in a r sin(o;t). For the disk flux, this yields 



1 + 



al 



sin(u; t) 



cos(2u; t) 



(1) 



The disk spectrum is approximately a blackbody spectrum at the temperature of the disk at 
the transition radius, and for the purpose of a simple analytical estimate, we assume that it is 
monochromatic at a disk photon energy Ej)(t) oc Tc(rt r ) oc (r tr £)~ 3//4 . For any observed photon 
energy E, we define the ratio 



e(i)^£/£ D (i)^ eo c: 3/4 . (2) 

Now, at any given time, a fraction f c ~ i (1 — j) of the disk radiation will intercept the 
quasi-spherical inner ADAF-like corona and serve as soft seed photons for Compton upscattering, 
producing the time-variable hard X-ray emission. Note that additional hard X-ray emission 
will be produced by Comptonization of synchrotron and bremsstrahlung photons produced in 
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the inner portions of the ADAF. However, since we do not assume significant changes of the 
structure of the inner ADAF in the course of the small-scale oscillations of the transition radius, 
this internal synchrotron and bremsstrahlung emission will constitute a quasi-DC flux component 
which does not contribute significantly to the variability properties of the source. Being mainly 
interested in the quasi-periodic variability, we do not consider this emission component here. The 
time-dependent Comptonization response function to the oscillating disk flux in this geometry can 
be parametrized as 



h(t, t ,€q) = 6(u) [au^ e -«/0 + B u^ 1 G ( 2R *W _ u ) } (3) 
(Bottcher & Liang 1998), where now t is the observing time, r is the time of soft photon emission 



at the transition radius, u = t — r, is the Heaviside function, A and B are normalization factors 
which are generally energy and time dependent. In the general case, the indices a and n, and 
the time "constant" (5 also depend on photon energy and time {(3 depends on time through the 
time-dependence of the extent of the corona and the disk photon energy). 

The observable energy flux at time t at photon energy E is then given by 



t oo 

F(E,t) = f d F D (E,t) + f c J drj dE D F d (E d ,t) h(t,r,e ), (4) 

-oo 

where fd = 1 — f c . For the purpose of the analytical estimate, we now focus on two regimes for 
the observed photon energies. At low energies, close to the peak energy of the disk emission, the 
observed light curve will be dominated by the direct disk emission (the first term on the r.h.s. of 
Eq. H|) and the contribution from reflection in the outer regions of the ADAF, proportional to 
B ~ c/(2i? tr [r]). For an ADAF-like density and temperature profile, we may set k = 1, i. e. the 
reflection light curve is flat over the light crossing time through the corona. The multiple-scattering 
term, proportional to A in Eq. ||, may be neglected at low photon energies. 

The observed signal at high photon energies, E ^> Ed will be dominated by the multiple- 
scattering term, i. e. we may set B = for high-energy photons. Since multiple Compton 
scattering results in a power-law spectrum, we parametrize the energy dependence of the 
normalization of this term as A = Ao e$ 7 x~ 37//4 , where 7 is the photon spectral index of the 
Comptonized spectrum, and x = 1 + sin(uT). The parameter a is generally energy dependent. 
However, we find that a constant value of a = 10 yields a reasonably good description of the 
multiple-scattering light curve for an ADAF-like temperature and density profile. To the same 
degree of approximation, we consider j3 as constant for a given energy channel, and note its scaling 
as (3 oc r tr . Furthermore, f3 increases with photon energy. Now, using the 5 function assumption 
for the disk spectrum in the integral in Eq. ||, the observed light curve reduces to 
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OO 

F(E,t)n f d F D (E,t)+f c F D>0 J du ^Aoe^x^u^e- 11 ^ + X - 2 B (e ) 6 - u^J } , (5) 

o 

where 5 = (37/4) + 1. 

Expanding all terms up to the first harmonic and writing 

F(E t) 

' — = Vo + Vi sm(ut) + rj2 cos(cjt) + 773 sin(2a; t) + 774 cos(2u; t) (6) 

Fd,o 

we find 



Vo = h (l + f ) + fc * ^ (l + \al) + /cA eo 7 (Vl>] + ^±l!a 2 [C 2 (a) + 5 2 (a)]) , (7) 

?7i = -/dflr - 2/ c — a r sin^ r - / c A e 7 5a r Ci(a), (8) 

m = 2fc—a r (1 - cos0 r ) + f c A (iQ y 5a r Si(a), (9) 

m = -hc-a 2 r sin 2 r - \f c A e^5{8 + 1) a 2 SC(a), (10) 

13 5 1 

774 = -^fda 2 . - -fc—a 2 r sm(f) r cos^ r + -/ c A)eo 7 5(£ + l)a 2 (S 2 [a] - C 2 [a]), (11) 

where 4> r = ^ r tr , and we have defined the integrals 

00 

Sl (a) . J * ,rt-/0sin W = frio) ~W» , (12) 





00 



S / A, „»-'e-»/"oos M = fr(«)=^, (13) 





00 



/ n f . «-i -u/fl ■ 21 \ P a T(a) f cos(aarctan[2/3w]) 1 

S 2 (a) - / d« - a *e "/^ J H = {l " (1 + [2 ^p )a/2 } ' ( 14 ) 

00 

„ , , /* , „_i _„/,g 2, \ f3°T(a) f cos(aarctan[2/3u;]) 1 

C 2 (a) = y d U U a l e «^cos 2 M = —^2 |l + , (15 ) 





00 



o^/ \ f j a-i -u/B ■ , \ /? a r(a) sin(aarctan[2/3w]) 

,SC(a) = Jduu a l e u ^ sm(uu) cos(^n) = ^ (1 + [2 /3a;] 2 )"/ 2 ( } 
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Defining the phases, relative to the Rt r oscillation, of the signal in a given energy band at the 
QPO frequency and the first harmonic, respectively, by 

F(t) = 770 + Ci sin(cjt + A Q ) + £ 2 sin(2w t + A H ) (17) 

with tanAQ = 772/771 and tanAg = 7/4/773, we find for a low-energy channel (with mean energy 
close to the peak energy of the disk emission spectrum): 

, . 2/ c ^ (1 - COS 4> r ) , . 

tan Aq(LE) « - n Jc " B \ - (18) 

, x fd + 1 /c- sin(2(/. r ) , 
tanA g LE 2 ^ J ■ (19) 

|/^(l-cos[20 r ]) 

Assuming that the direct disk emission (oc /^) is strongly dominating at this energy, and that both 
773 and 774 are negative in this case, we have Aq(LE) ~ 0, and Ag(LE) ~ —ir/2. 

For high photon energies, dominated by multiple Compton scattering, we find 



Aq(HE) = -aarctan(/5u;) + k 7r, (20) 

Ag(HE) = a arctan(2 /3cj) + k\iT, (21) 

where ki S (0, —1, 1) are determined by the signs of 771, . . . , 774. With the definition of the phases 
A in Eq. [l7], the phase lags are given by A<f> = A(LE) — A (HE). Throughout this paper we adopt 
the convention that positive phase and time lags correspond to hard photons lagging the soft ones. 

Intuitively, phase and time lags between the disk-radiation dominated low-energy photons and 
the Compton-upscattering dominated high-energy photons are determined by the ratio between 
the QPO period Tqpo and the time required for soft photons to reach the inner ADAF region 
and be Compton upscattered to hard X-ray energies. If this light-travel and diffusion time is 
longer than half a QPO period, but less than the QPO period itself, the phase lag at the QPO 
fundamental frequency — which physically is still a hard lag — appears as a soft lag due to the 
periodicity of the light curves. The same argument holds for the n th harmonic of the QPO, where 
the period now corresponds to l/(n + 1) of the QPO fundamental period. 

For a very large transition radius, the direct disk emission may not contribute significantly to 
the X-ray flux. In that case, the phase lag between two energy channels will be dominated by the 
difference in diffusion time scales (3, and we expect 

A0QPO ~ a(arctan[/?ju/] — arctan[/3jw]), (22) 
harm. ~ a(arctan [2/3jO>] — arctan[2/3jw]). (23) 
for two energy channels E{ < Ej. 
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3. Monte Carlo simulations 

The estimates derived in the previous section were based on several major simplifications in 
order to keep the problem analytically trackable. In order to test the validity of our approximations 
and to provide results under more realistic assumptions, we have simulated the radiation transfer in 
our oscillating ADAF/disk model system, using our time-dependent Monte-Carlo Comptonization 
code. For a detailed description of the code and its capabilities, see Bottcher & Liang (1998, 1999). 

In our simulations, we approximate the accretion disk emission by a blackbody spectrum 
with the temperature of the disk inner edge at any given time, and its time-dependent luminosity 
determined as described in the previous section. A fraction f c of the time dependent disk emission 
enters the spherical Comptonizing region at its outer boundary and serves as seed photon field for 
Comptonization. The ADAF is characterized by an r~ 3 / 2 density and r~ l temperature structure, 
and the electron temperature is normalized to kT e (R s ) = 500 keV, i. e. the electrons become 
relativistic at the event horizon. The event horizon is treated as an absorbing inner boundary. 
The corona has the specified radial Thomson depth of r^ DAF when the ADAF/disk transition 
is located at r tr . As r tr oscillates, we leave the density and temperature structure of the ADAF 
inside Rt T (t) unchanged and set the coronal electron density equal to zero outside Rtv(t). 

The resulting light curves, generally consisting of the sum of direct disk emission and the 
Comptonized emission from the corona, are sampled in 5 photon energy bins and over 512 time 
steps of At ^ 0.05 Tqpo- For several test cases, we have run identical problem simulations with 
different time steps, in order to verify that our results are independent of the time step chosen. 



The energy-dependent light curve are then Fourier transformed, using an FFT algorithm (Press et 
al. 1992| ), and the power density spectra and phase lags are calculated. 



In order to test our Monte Carlo code against the analytical estimates derived in the previous 
section, we did a series of simulations in which the accretion disk temperature at the equilibrium 
transition radius and the QPO frequency were artificially chosen artificially constant among 
different simulations, and at a value such that the emission in the 2-5 keV reference channel was 
always dominated by direct disk emission, while the highest two energy channels (15 - 40 keV 
and 40 - 100 keV, respectively) were dominated by Compton-upscattered radiation from the 
corona. In that case, Eqs. (18) - (^) with constant values of a and /3/ r tr may be used in order to 



approximate the expected phase lags. In Fig. ||, the phase lags at the QPO fundamental and first 
harmonic frequencies as measured in our simulations are compared to these analytic estimates. 
The figure demonstrates that the two approaches are in excellent agreement, indicating that our 
numerical procedure reliably reproduces the time-dependent radiation transport properties of the 
model system. 
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4. The 0.5 10 Hz QPOs in GRS 1915+105 

As briefly mentioned in the introduction, the variable-frequency QPO at 0.5 - 10 Hz observed 
in the low-hard state of GRS 1915+105 exhibits a very peculiar phase lag behavior ( |Reig et~al 



2000] |Lin et al. 20"00| ): If the source has a rather hard (photon index 7 ^ 2.7) photon spectrum 



and low X-ray flux, the QPO frequency is low /qpo ^ 2 Hz, and the phase lags associated with 
both the QPO fundamental and the first harmonic frequencies are positive. As the X-ray flux 
increases and the photon spectrum becomes softer, the QPO centroid frequency increases, and the 
phase lag at the QPO fundamental frequency decreases and changes sign at /qpo ~ 2.5 Hz. At 
the same time, the phase lag associated with the first harmonic frequency increases. As the source 
becomes more X-ray luminous and the spectrum becomes softer, the QPO frequency increases up 
to ~ 10 Hz, and the phase lag associated with the QPO fundamental continues to decrease, while 
the phase lag at the first harmonic frequency remains positive and does not show any obvious 
correlation with spectral parameters or the QPO frequency. Spectral fits to GRS 1915+105 with 
a disk blackbody + power-law model to different spectra states along this sequence revealed that 
the inner disk temperature increases from ~ 0.7 keV to ~ 1.5 keV as the QPO frequency increases 
from 0.5 to 10 Hz flMuno et al. 1999|) . 



In order to model this behavior, we parametrize the disk temperature at the transition 

-3/4 



tr 



radius, TD(r tr ), and the QPO frequency as power-laws in r tr . Assuming TD(r tr ) oc r 
a simple power-law relation spanning the observed ranges of 0.7keV < kT-Q < 1.5keV and 
0.5 Hz < /qpo < 10 Hz ( Muno et al. 199S| ) yields a scaling /qpo oc T^j oc r t ~ 3 . This scaling may 



indicate that the transition radius oscillations are related to a modulation of the disk evaporation 
process, responsible for the transition to the inner ADAF, by a secular instability. The frequency 
of such modulations might be proportional to the inverse of the evaporation time scale, 7" ev ap> 
which is expected to be proportional to the disk surface flux, dL/dA oc T^j oc r^ 3 . Thus, the 
assumed scaling laws of the disk temperature and QPO frequency with the transition radius are 
physically plausible. 

For the high-QPO-frequency end of the sequence mentioned above, we choose a transition 
radius of r tr = 6 • 10 8 cm, corresponding to ~ 700 R s for a 3 M@ black hole or the theoretical limit 
of 340 R s found by Meyer et al. (2000) for a ~ 6M Q black hole (unfortunately, the mass of the 
central object in GRS 1915+105 is not known due to the lack of an optical counterpart). The 
QPO frequency (i.e. the transition-radius oscillation frequency) is 10 Hz, and the disk temperature 
at the transition radius is 1.5 keV. Assuming M = 6M and rt r = 340 R s , this is in agreement 
with the cool, gas-pressure dominated disk model for a -1 / 5 (L/0.057Le ( j ( j) 3 ^ 10 ~ 1 5 where a is the 
viscosity parameter. 

Using the Monte-Carlo simulations described in the previous section, we calculate the photon- 
energy dependent light curves, Fourier transform them, and calculate the Fourier-frequency 
dependent phase lags from the cross-correlation functions. For comparison with the results of Lin 
et al. (2000) we sample the light curves in the energy channels 0.1 - 3.3 keV [1], 3.3 - 5.8 keV 
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[2], 5.8 - 13 keV [3], 13 - 41 keV [4], and 41 - 100 keV [5], and focus on the phase lags between 
channels [4] and [2]. We perform a series of 6 simulations, each with different values of r tr , implying 
different QPO frequencies and disk temperatures according to the scaling laws quoted above. 

Two representative examples of the resulting power and phase lag spectra are shown in Figs. 
H and |3] for a high-frequency (6.7 Hz) case with alternating phase lags and a low-frequency (1.1 Hz) 
case with positive phase lags at both fundamental and first harmonic frequencies, respectively. 

The results of the complete series of simulations are listed in Table [l] and illustrated in Fig. 
||. From Fig. || we see that this model reproduces the key features of the peculiar QPO phase lag 
behavior observed in GRS 1915+105, in particular the change of sign of the phase lag associated 
with the QPO as the QPO centroid frequency increases, and the alternating phase lags between 
QPO fundamental and first harmonic at high QPO frequencies. We point out that the measured 
values of the hard lags at the first harmonic frequency for /qpo 2 Hz are significantly smaller 
than predicted by our model. However, this is expected because, for realistic signal-to-noise ratios, 
the (random-phase) Poisson noise contaminating the actual data always suppresses high phase lag 
values, but leaves small values of \A(j>\ almost unaffected (Zhang 2000, private communication). 

Table [l] and Fig. |I] also show the predicted phase lags at the second harmonic of the QPO, 
which has not been observed yet. We expect that, if a second harmonic is found in any future or 
archived observation of GRS 1915+105, it should be associated with a negative phase lag, at least 
for /qpo £ 2 Hz. 

We need to point out that the more complicated case discussed in this section is not directly 
comparable to the test cases illustrated in Fig. [I] since we are now varying the accretion disk 
temperature and the QPO frequency when varying the transition radius. This also changes the 
parameters a and f3 entering the analytical estimates as a function of r tr in a non-trivial way. 

5. Summary and conclusions 

We have investigated the time-dependent X-ray emission from an oscillating two-phase 
accretion flow, consisting of an outer, cool, optically thin accretion disk, and an inner ADAF. 
Based on this model, we are proposing an explanation for the peculiar phase lag behavior associated 
with the variable-frequency 0.5 - 10 Hz QPOs observed in GRS 1915+105. In particular, the 
changing sign of the phase lag at the QPO fundamental frequency as the QPO frequency increases, 
and the alternating phase lags between QPO fundamental and first harmonic frequencies in the 
case of a high QPO frequency are well reproduced by this model. The relation between QPO 
frequency, transition radius and inner disk radius can be interpreted physically, if the QPO is 
triggered by a modulation of the transition radius by a secular instability affecting the disk 
evaporation responsible for the disk/ADAF transition. Based on our results for GRS 1915+105, 
we predict that, if a second harmonic to the variable-frequency QPO at /qpo ^ 2 Hz is detected, 
it should be associated with a negative phase lag. 
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We found that oscillations at the transition radius can naturally lead to apparent soft lags 
associated with the QPO frequency, if the photon diffusion time related to the production of hard 
X-rays is longer than half a QPO period. We point out that it is very unlikely that the same 
mechanism is responsible for the alternating phase lags observed in the 67 mHz QPOs and its 
harmonics since this would require a transition radius at rtr ^ 10 5 R s . This large size of the ADAF 
is implausible for a non-quiescent state. The 67 mHz QPO alternating phase lag pattern has been 
observed in a high- flux state of GRS 1915+105, in which we would generally expect the disk to 
extend down to rather small radii. 

One caveat of the model we have proposed here is that it requires a rather large accretion 
rate in order to produce the observed disk luminosity from GRS 1915+105. This, in turn, requires 
either that a large fraction of the accreted mass is lost to the collimated outflow of the radio 
jets, or that the ADAF has a very low radiative efficiency in order not to overproduce the hard 
power-law emission from the hot ADAF. However, since GRS 1915+105 is known to have very 
powerful, superluminal radio jets ( Mirabel fe Rodriguez 1994 ), it is plausible to assume that, 



indeed, a significant fraction of the mass accreted through the outer disk is not advected through 
the ADAF, but is powering the radio jets. 

The work of MB is supported by NASA through Chandra Postdoctoral Fellowship Award 
Number PF 9-10007, issued by the Chandra X-ray Center, which is operated by the Smithsonian 
Astrophysical Observatory for and on behalf of NASA under contract NAS 8-39073. 
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Table 1. Results of our Monte-Carlo simulations to reproduce the phase lag 
behavior of GRS 1915+105. The phase lags Acfi [rad] are the hard lags between 
channels [4] (13 - 41 keV) and [2] (3.3 - 5.8 keV). 



r tr [cm] 


/qpo [Hz] 


kT D [keV] 


A^qpo 


A^l.harm. 


harm. 


1.7 x 10 9 
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+0.11 


+0.13 


-0.58 


1.3 x 10 9 


1.1 
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+0.11 


+0.20 


+0.20 


1.0 x 10 9 
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1.0 


+0.06 


+1.71 


-3.13 


8.0 x 10 8 
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1.2 
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+2.54 


-2.57 


7.0 x 10 8 


6.7 


1.4 


-0.13 


+2.64 


-2.59 


6.0 x 10 8 


10.0 


1.5 


-0.17 


+2.54 


-2.92 
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Fig. 1. — Comparison of the analytical estimates for the phase lags at the QPO fundamental and 1. 
harmonic frequencies to the results of numerical Monte-Carlo simulations. Solid curve: analytical 
approximation for the lag at the QPO fundamental frequency; dot-dashed curve: approximation for 
the lag at the 1. harmonic; filled symbols: phase lags at the QPO fundamental frequency resulting 
from simulations; open symbols: phase lags at the 1. harmonic from simulations. Parameters: 
/ QPO = 10 Hz, t4 DAF = 2, kT D (r tT ) = 0.5 keV. The energy channels are [1] 0.1 - 2.0 keV; [2] 2.0 
- 5.4 keV; [3] 5.4 - 15 keV; [4] 15 - 40 keV; [5] 40 - 100 keV. 
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10 

f QPO [Hz] 

Fig. 2. — Power spectra (top panel) and phase lags with respect to the 3.3 - 5.8 keV channel (bottom 
panel) of our model simulations at a high QPO frequency of /qpo = 6.7 Hz, with r tr = 7 x 10 8 cm. 
Disk temperature: kTj^(r tT ) = 1.4 keV. Dashed vertical lines indicate the location of the QPO and 
its harmonics. The phase lag at the QPO fundamental frequency is negative, while the lag at the 
first harmonic frequency is positive. For the second harmonic, a negative lag is predicted. 
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Fig. 3. — Power spectra (top panel) and phase lags with respect to the 3.3 - 5.8 keV channel (bottom 
panel) of our model simulations at a low QPO frequency of /qpo = 11 Hz, with r tr = 1.3 x 10 9 cm. 
Disk temperature: kT-£){r tT ) = 0.8 keV. Dashed vertical lines indicate the location of the QPO and 
its harmonics. The phase lags at both the QPO fundamental and the first harmonic frequencies 
are positive for the 13 - 41 keV channel. This simulation predicts a positive phase lag also at the 
second harmonic. 



-16- 



T3 
03 



0.3 

0.2 

0.1 

0.0 

-0.1 

-0.2 

-0.3 
1.6 

1.3 

1 

0.7 



E 15 
o 

00 

° 10 



> 

CD 







♦ ♦• . 



■■"■•»<... 



♦ A ^QPO 

-a 0.1 x Harmonic 
-•0.1 xA^ Harmonic 



\ 



8 



10 



f Q Po [Hz] 

41 keV and the 3.3 



Fig. 4. — Phase lags between the 13 - 41 keV and the 3.3 - 5.8 keV energy channel at the 
QPO fundamental and 1. harmonic frequencies (upper panel) resulting from our Monte-Carlo 
simulations, with inner disk temperature as a function of QPO frequency (middle panel) determined 
from observations of GRS 1915+105 and disk/ADAF transition radius (lower panel) determined 
through To(rtr) oc r tr 3 ^ 4 . The simulations properly reproduce the changing sign of the phase lag 
as the QPO frequency increases, and alternating phase lags for large QPO frequencies, as observed 
in GRS 1915+105. 



